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The Giessen Bohzmann-Uehhng-Uhlenbeck transport model is extended 
and applied to the p-nucleus interactions in a wide beam momentum range. 
The model calculations are compared with the experimental data on p- 
absorption cross sections on nuclei with an emphasis on extraction of the 
real part of an antiproton optical potential. The possibility of the cold 
compression of a nucleus by an antiproton in-flight is also considered. 

PACS numbers: 25.43. -|-t;21.30.Fe;24.10.Lx;24.10.-i 

1. Introduction 

Antiproton interactions with nuclei are of big interest since they deliver 
the information on in-medium pN interactions [I] related to the antiproton 
optical potential. The real and imaginary parts of p optical potential close 
to the nuclear centre are still poorly known ^ |3]. As shown in [1]-[Z]) a 
deep enough real part, Re(V^pt) = —(150 — 200) MeV, might cause sizeable 
compressional effects in p-doped nuclei. 

Here, we present some results of the microscopic transport simulations 
within the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model [8j. Our 
goal is twofold: (i) to determine the p optical potential by comparison with 
the data on p absorption cross sections on nuclei; (ii) to evaluate the proba- 
bility of p annihilation in a compressed nuclear zone for the beam momenta 
0.1-10 GeV/c. 

The model is described in Sec. II. In Sees. Ill and IV, the calculations 
of p absorption cross section on nuclei and of the dynamical compression 
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of the O nucleus by moving antiproton are presented. The results are 
summarized in Sec. V. 



2. The GiBUU model 

The GiBUU model [8] solves the coupled set of kinetic equations for 
different hadrons {N, N, A, A, tt...). These equations describe the time 
evolution of a system governed by the two-body collisions and resonance 
decays. The Pauli blocking for the nucleons in the scattering final states 
is accounted for. Between the two-body collisions, the particles propagate 
according to the Hamiltonian-like equations in the mean field potentials 
defined by the (anti)baryon densities and currents. 

We use the GiBUU model in th e relativist ic mean field mode. The 
single-particle energies = + ^Jp^'+m*^ of the nucleons {i = N) 

and antinucleons {i = N) depend on the effective mass m* = rrii + Si, 
where Si = g^ia is the scalar field, kinetic three-momentum p* = p — Vj, 
and vector field V/^ = g^jii^'^ + QpiT^p^^ + eiA^^. Here, a, w, p and A 
are, respectively, the isoscalar-scalar, isoscalar-vector and isovector-vector 
meson fields, and the electromagnetic field. 

The meson- nucleon coupling constants gaN-, g^N and gpN are taken 
from the NL3 set of parameters of a non-linear Walecka model. The 
meson-antinucleon coupling constants are motivated by the G-parity trans- 
formation, however, allowing for their strengths to be rescaled by a factor 
< < 1 as = igaN, 9ujN = -C9loN and = S,gpi^. 

The following two-body collision processes involving an antinucleon are 
implemented in the model: elastic scattering and charge exchange NN 
NN, inelastic production NN BB + mesons, and annihilation NN 
mesons. Further details of the model can be found in [51 El [TDl [TH [T2] and 
in refs. therein. 



3. Antiproton absorption and annihilation on nuclei 

Fig.[Tjshows the antiproton absorption cross section on ^^C. The GiBUU 
calculation without nuclear part of the antibaryon mean field (f = 0) is in 
a very close agreement with the Glauber model prediction [13ji. 

The attractive real part of an antiproton optical potential 

Re(Kpt) = S, + V^+ ^ ^ (1) 



^ Pure G-parity transformed nuclear fields are obtained with ^ = 1. But this leads to 
an unrealistically deep real part of the p optical potential, Re(Vopt) — —660 MeV. 

^ Except for momenta less than ~ 0.2 GeV/c, where the Coulomb potential causes 
the deviation. 
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Fig. 1. (Color online) Beam momentum dependence of the antiproton absorption 
cross section on ^^C. GiBUU calculations are represented by the curves denoted 
by the value of the scaling factor ^ of meson-antibaryon coupling constants. The 
results obtained using the Glauber model [13] are shown by thin solid line. The 
calculation with ^ = without annihilation is shown by dotted line. Experimental 
data are from 

bends the p trajectory towards the nuclear centre and increases the p ab- 
sorption cross section. One can see from Fig. [U that the finite values of 
the scaling factor ^ ~ 0.2 are needed to describe the antiproton absorption 
cross section data at piab = 470 — 880 MeV/c measured at KEK [2]. The 
best fit of the data [2j on p absorption cross sections on the "^^C, ^'''Al and 
^■^Cu nuclei by GiBUU calculations is reached with ^ ~ 0.21 it 0.03, which 
produces Re(yopt) — —(150 ± 30) MeV in the nuclear centre. 

The imaginary part of the p optical potential can be calculated as 



where VpN is the relative velocity of the antiproton and a nucleon; cr^f is 
the total pN cross section including the in-medium effect of Pauli blocking 
for the final state nucleon in the NN — )• NN channel; /Oat is the local 
nucleon density. The averaging in Eq.([2]) is taken with respect to the nucleon 
Fermi distribution. Using ([2]), we obtain Im(V^pt) — —(105 it 5) MeV in 
the nuclear centre, where a small uncertainty is due to different considered 
nuclei. More details on the extraction of an antiproton optical potential by 
GiBUU calculations can be found in |12] . 



Im(Kpt) = "2 < "^pW^Sf > PN 



(2) 
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Fig. 2. (Color online) Nucleon (thick solid lines) and antiproton (thin solid lines) 
densities along the z-axis passing through the nuclear centre (x, y, z) = (0, 0, 0) for 
the p^^O system at different times. The antiproton has been initialized at the coor- 
dinates (0, 0, —2.5) fm with momentum components (0, 0, 0.3) GeV/c. The scaling 
factor of the antiproton-meson coupling constants ^ = 0.22 is used. Annihilation 
is switched-off. 



4. Dynamical compression of a nucleus by an antiproton 

Fig. [2] shows the time evolution of nucleon and antiproton densities for 
the case of an antiproton initialized at the nuclear periphery with momen- 
tum of 0.3 GeV/c directed towards the nuclear centre. The calculation has 
been done without the NN annihilation channel in the collision term, but 
allowing for the NN and A^A^ scattering. We observe that the antiproton 
attracts surrounding nucleons catching them into a potential well of about 
70 — 100 MeV depth. The nucleon density bump moves with the p or slightly 
behind it due to some delay in the reaction of the nucleon density on the 
perturbation created by the antiproton. 

The compression process of Fig. [2] has to be complemented with the 
calculation of an antiproton survival probability 

Psurv(t) =exp|- J dt'rann(0| • (3) 

where Fann = Pn {vpNO' ann) is the autiproton annihilation width and Cann is 
the pN annihilation cross section (other notations are the same as in Eq. ([2]) ) . 
In the case of the process shown in Fig. [21 the antiproton survives with the 
probability Psurv ^ 10~^ at the time 10 fm/c when the nucleon density max- 
imum is reached. The experimental datection of the nuclear compression by 
the antiproton is only possible, if p annihilates in the compressed zone of a 
nucleus. The corresponding probability is Pcompr(Pc) = Psmv{ti) — Psurv{i2) ■, 
where ti < t2 are the times delimiting the time interval, when the maximum 
density of a nuclear system exceeds some preselected value pc- In concrete 
calculations, we set pc = 2po with po = 0.148 fm~^ being the nuclear satu- 
ration density. 
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Fig. 3. The probability of antiproton annihilation at the nueleon density exceeding 
2/9o as a function of the beam momentum for p^^O collisions. 

In the case of a real antiproton-nucleus collision, we apply, first, the 
GiBUU model in the standard parallel ensemble mode to determine the 
coordinates r and momentum p of the antiproton at its annihilation time 
moment event-by-event. Next, we initialize p at (r, p) and run GiBUU with- 
out annihilation. This allows us to compute i-*compr for a given annihilation 
event, i.e. to determine the probability that this event will take place in a 
compressed nuclear zone. Finally, we calculate the cross section of a p an- 
nihilation in the compressed zone by weighting with the impact parameter 
as (Jcompr = /o°° db 2Trb PcompriPc b), where the upper line denotes averaging 
over annihilation events with a given impact parameter b. 

In Fig. [3l we present the probability of p annihilation in a compressed 
zone given by the ratio o"compr/o"ann; where cJann is the total annihilation 
cross section of p on the oxygen nucleus. The rise of the ratio with the 
beam momentum between about 1 and 3 GeV/c is caused by opening the 
pion production channels NN — )• NNtt^tt...), which leads to more inten- 
sive stopping of an antibaryon before annihilation, and, therefore, to larger 
compression probabilities. 

5. Summary and conclusions 

We have performed the microscopic transport GiBUU simulations of p- 
nucleus collisions. The antiproton mean field potential has been described 
applying the non-linear Walecka model supplemented by appropriate scaling 
of the meson-antinucleon coupling constants. 

The extracted antiproton optical potential from comparison with the 
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KEK data on p absorption cross section below 1 GeV/c is Vopt = —(150 it 
30) — i(105ib5) MeV in the nuclear centre, which is comparable with the well 
known phenomenological values [21 [3]. However, the BNL and Serpukhov 
data on p absorption above 1 GeV/c require much deeper real part of about 
—660 MeV, close to the G-parity value. It is important, therefore, to perform 
the new measurements of the p absorption cross section on nuclei above 1 
GeV/c at FAIR. 

The probability of the antiproton annihilation in a compressed nuclear 
zone is maximal at the beam momenta below 0.5 GeV/c. However, it is 
also big enough (~ 10^^) at the beam momenta from 3 to 10 GeV/c. The 
possibility of additional triggering, e.g. on a high-momentum proton, would 
prefer this beam momentum range with respect to the low beam momenta 
for the study of annihilation events in the compressed nuclear zone. 
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